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Abstract 



Top quarks produced in pairs are predicted to experience spin correlations. Due to the 
large tt statistics expected for the LHC, it should be possible to search for new physics 
effects in angular variables sensitive to these correlations. We investigate, for a gen- 
eral two-Higgs Doublet Model (2HDM), the charged Higgs boson decay of the top quark 
through the channel t — ► bH + — ► br + v r . Analytic results are presented on the spin an- 
alyzing coefficients for this decay mode. We then explore in some detail the correlation 
phenomenology in the Type II 2HDM. Finally we present a hadron-level Monte Carlo 
analysis, illustrating distributions in azimuthal angles which are sensitive to correlations 
in the transverse plane. These observables are accessible also in the r channel, and are 
therefore particularly interesting for analyzing the t — * bH + decay. 



1 Spin correlations in tt hadroproduction 

When top quarks are pair-produced in a hadronic collision, the spin projections of the t and 
the t can be correlated to a certain degree when a suitable basis is chosen for spin quantization. 
There exists an extensive earlier literature on this topic; see for example [1-3], or the more 
recent review [4]. 

The degree of tt spin correlation can be expressed as 



where a, b are interpreted as spin quantization axes for on-shell t, t. To determine the value of 
this observable, as a function of the invariant mass of the top pair m t j, one folds the parton- 
level correlations Cy of each partonic subprocess with the parton distribution functions (pdfs). 
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Here we work exclusively in the helicity basis, defined by a = p t and b = pj, for which the 
simple threshold expressions C qq - = — 1/3 and C gg = 1 hold. In the UV limit CV,- — > — 1 due to 
helicity conservation. From an NLO QCD calculation [5], a total helicity correlation C = 0.326 
has been determined for the LHC (y/s = 14 TeV). At leading order the corresponding value is 
C = 0.319. These numbers were obtained using m t = 175 GeV and CTEQ6.1 pdfs. 

As a special case of the correlation given by eq. ([Tj), it is possible to consider the summed 
contributions from correlations along three orthogonal axes. This corresponds to 



At the LHC the expected value for this observable is V = -0.219 (-0.212) at NLO (LO) [5]. 
The ample ti statistics expected with cr^ LO ~ 900 pb at y/s = 14 TeV will allow the LHC 
experiments to measure these correlations for the first time. We note in passing that, for 
these energies, the degree of correlation in the helicity basis can be substantially increased by 
sacrificing some statistics through a cut on the maximal m t j [2]. 

2 Polarized top quark decay and angular correlations 

In the rest frame of a decaying (on-shell) top quark, with the spin directed along the z-axis, 
the differential decay rate can be written [6] 



with cos^j = pi ■ z for each decay product i with momentum direction pj. Each decay product 
also has a spin analyzing coefficient ai, which follows from the Lorentz structure of the Wtb 
coupling. The ai are obtained by integrating the polarized decay matrix element over the 
appropriate phase space. Performing this calculation, one obtains for the SM with m t = 172.6 
GeV that aw = — &b = 0.39, ai = a^ = 1, and a n = a u = —0.34. The analytic expressions are 
summarized in table [TJ 

In order to utilize the correlations described in sect. [TJ through decay distributions, one 
decay product from each top quark are combined into a doubly differential distribution 



Terms linear in the cosines are forbidden by parity invariance. The angles 0i(6j) are determined 
in the respective rest frames of t(t), as the angles to the directions of motion of t(i) in the 
overall ti CM frame. In analogy with eq. ([2]), a similar angular distribution in cos 6^- = pj • pj 
is obtained for the P-type correlation 
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Table 1: Spin analyzing coefficients ctj for different decay products within the SM and 
for a charged Higgs decay of the top quark (t — > bW + /H + — > bl + b>i), or equivalently 

(t -» -» bdu). 

Analyzing Coefficient a% for decay channel 

particle W + (u = mly/m 2 ) H + (£ = m 2 H+ /m 2 ) 
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As noted in [7], the distribution given by eq. (jSJ) is less sensitive to acceptance loss by phase- 
space cuts than that given by eq. (BJ. 

It is possible to measure this type of distributions experimentally to obtain information on 
C, T> (assuming SM decay) or cti, a.j (when the correlation can be determined in an independent 
manner). Correlation measurements at the LHC were studied previously for ATLAS [7], with 
an expected ~few percent precision on determining C and T> in the dilepton channel («/«/ = 1) 
using 10 fb _1 of data. Such an analysis would not be statistics limited, inviting us to study 
rare processes related to new physics together with top quark spin correlations. This will be 
the topic of the next section. 



3 Top quark decay to charged Higgs bosons 

Should the top quark have non-standard couplings, changes can occur either to the production- 
level correlation, and/or to the decay distributions. A general account of anomalous Wtb 
couplings in single top production is given in [8]. If they exist, the same couplings are of course 
relevant to top decay. 

We shall take a different route and discuss instead a charged Higgs decay of the top quark. 
This mode could become important in two-Higgs Doublet Models (2HDM) when < m t — 
rrib- From LEP, the direct limit on the mass of a charged Higgs boson is > 79.3 GeV at 
95% CL [9], assuming only the decays H + — ► t + v t and H + — > cs are possible. Parametrizing 
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Figure 1: Branching fraction for top quark decay t — ► bH + in the 2HDM (II) for different values 
of rriH+ ■ 

the H + interactions with fermions as 



we determine the ai for the decay t — > bH + — > bl + ui as a function of the couplings A, B, 
and C. The calculation is performed in the narrow- width approximation for on-shell t(t), 
assuming independent decay channels. Further details are described in [10]. The results are 
given in table [TJ For the decay through H + , the coefficients are seen to depend on two universal 
factors. Since f(£,A,B) is very close to unity except in the limit m#+ — > m t , when anyway 
BR (t -> bH + ) -> 0, it will simply be set to 1. The coupling factor (A 2 - B 2 )/(A 2 + B 2 ) is 
important, and will be further discussed below. Our results agree with those in [11], where also 
0(a s ) corrections to an are calculated. 

Motivated by minimal supersymmetry, the most widely studied 2HDM is the Type II, where 
one doublet is coupled exclusively to up-type fermions, and the other to down-type fermions. 
In this model the couplings have values A = m u cot/3, B = m^tan/5, and C = mi tan. P, where 
tan/3 is the ratio of the two doublets' vacuum expectation values. Since the Higgs couples 
proportionally to the fermion masses, we are concerned only with couplings to third generation 
particles in the following. QCD corrections to the Yukawa couplings are taken into account by 
renormalizing the fermion masses in the MS scheme at the scale m#+ [12]. Fig. [T] shows the 
resulting branching ratio t — > bH + as a function of tan /3 for different . 

In fig. [2] we illustrate the tan/3 dependence of cn& for the 2HDM (II). As seen in the figure, 
the efficiency to analyze the top quark spin is highest in the limits of small and large tan/3, 
where it can reach unity. The b quark (or the H + itself, since = — a&) is therefore the 
most powerful analyzer in the 2HDM. This is a difference compared to the SM, where it is 
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Figure 2: Spin analyzing coefficient a b for top decay in the 2HDM (II) and the SM. Fermion 
masses in the Higgs couplings are evaluated at the scale m H + = 100 GeV. The solid curve 
shows the pure 2HDM (II) result, while dashed curves correspond to scenarios with large SUSY 
corrections e& = — e' t = ±0.01. 

the charged lepton (or, equivalently, the down-type quark from W — > ud) which carries unit 
efficiency 

We consider also SUSY QCD corrections to the H + couplings, which occur in the 2HDM 
(II) of the MSSM. The tan /^-enhanced corrections to the down-type Higgs doublet coupling are 
studied in [13, 14]. They can be correctly resummed and evaluated to all orders by replacing 

tan/? 

m b t&nfj -> m b — 

1 + e& tan p 

m t cot (3 — > m t cot — e' t tan 0) 

in the couplings. Here €b,e' t are functions which fulfill |e,tan/3| < 1. As can be seen in fig. [21 
the effects on the oij are small; in the transition region with intermediate tan /?, they amount 
to at most a 10 — 20% correction. The effect cancels completely in the high tan/? limit. We 
note that similar SUSY corrections apply also in the limit of large cot (3. The conclusions on 
OLi should be the same for this case. 

To compare the angular correlations in the SM with those in the 2HDM (II) we display 
in fig. [3] analytic results on the D-type distributions from eq. ((SJ). The results are given for 
the process pp —>■ ti — > bbT + u T ud, where rv T in the case of 2HDM originates from the decay 
H + — ► t + v t which dominates over a broad tan f3 interval. The two combinations (r, d) (optimal 
for SM) and (b, d) (optimal for 2HDM) are chosen to illustrate the maximal differences. In the 
figure the two extreme cases for tan/? are shown with the SM expectation, and with a flat 
distribution corresponding to no ti spin correlations. 
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Figure 3: Truth distributions for different final state particles The lines correspond to 

the SM with (solid black) and without (black dotted) ti spin correlations, the 2HDM (II) for 
tan/5 = 1 (dashed blue), and for tan/3 = 50 (dotted red). The results were evaluated using 
m H+ = 80 GeV. 

4 Monte Carlo Simulations 

There are certain limitations in the treatment presented in the previous section, namely that the 
t decay of the H ± prevents reconstruction of the ti rest frame since there are too many neutrinos 
in the final state. This reconstruction is necessary in order to construct the angular distributions 
discussed above. To address this issue, we consider correlations in the plane transverse to the 
beam [10]. The angular variable analogous to cos% in this case is p_u ■ p±j = cos Atfiij, where 
P±i;P_i_i are the transverse direction vectors in the transverse rest frames of t(i). Introducing 
the notation A0j for the azimuthal angle measured to a fixed axis, for particle i in the transverse 
rest frame, A0y = A0j — A<pj. Using these variables, we expect a distribution of the form 



dN 1 



iVdcos(A0 f - A(j)j) 2 



1 + VonOLj cos(A0j - A(f)j) . (7) 



The transverse correlation at the LHC we determine numerically to have the value T>' = —0.193 
at LO. Comparing figs. [3] and HI their similarity shows that most of the information on the ti 
spin correlations is indeed accessible in the transverse plane. 

To test if these parton level results apply also at the hadron level we perform a Monte Carlo 
study of the signal pp — ► ti — > bb(H ± — » t ± u t )(W t — ► qq), with r — > hadrons. In addition we 
generate the same final state for the SM process. The matrix element for the 2 — > 6 process is 
generated with MadEvent [15] to preserve the spin correlations. PYTHIA [16] is then applied 
for parton showering, hadronization and generation of the underlying event. The hadronic r 
decay is treated by TAUOLA [17]. 

A simple spin correlation analysis is employed as follows: We require all particles to have 
\r)\ < 5 to account for the detector region. An exclusive mode k± jet finding algorithm is used, 
with the separator d cut = 20 GeV. Jets with \t]\ < 2.5, which are matched to a true b quark or 
r lepton, are identified and flavour "tagged". The identification of exactly one such r jet and 
two b jets is required. W and t candidates are constructed from the jj and jjb combinations 
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Figure 4: Azimuthal distributions for different final state particles (i, j). Color coding as in 
fig. 3. All results for m H + = 80 GeV. 
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Figure 5: Truth distributions for different final state particles Color coding as in fig. 3. 

All results were obtained for m H + = 80 GeV. 



which minimizes the \ 2 difference to the known W and t masses. It is required that the W 
(t) candidates are within 10 (15) GeV of the true mass values for the events to be retained for 
further analysis. 

Since one of the final state b jets is assigned to the t decaying in the fully hadronic mode 
(SM), the other b jet can be unambiguously assigned to the prospective H + side of the event 
(where also the r originates from). We reconstruct the transverse momentum of the t on this 
side using the combination p f ± = pj_ ct + p^ et + p™ lss . 

Unfortunately, neither of the most sensitive distributions, presented in fig. HI are experimen- 
tally accessible. This is because there is no lepton originating from the SM side of the event, and 
the purity in identifying the down-type quark from W — > qq' is poor. Hence the less efficient 
b jet (|a^ M | = 0.39) must be used. In fig. Owe show the analytic form of the two accessible 
distributions (r, b) and (b, b) . From the MC events we can evaluate the same observables at the 
hadron level. The result is given in fig. [6] for = 80 GeV. The choice of using this very low 
value for m#+ is simply to isolate the genuine spin effects from the kinematic differences. 
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Figure 6: Jet level distributions for different final state jets Colour coding as in fig. 3. 

Results obtained for m H + = 80 GeV. All distributions are normalized individually. 

Since we choose to work with normalized distributions, the different histograms in fig. O are 
not normalized to any particular branching ratio t — > bH + . The bin- by-bin fluctuations can be 
taken as a measure of the statistical variations in the number of events corresponding roughly 
to fO fb -1 , remembering the crude analysis employed. That the azimuthal distributions have 
the same advantage as the P-type correlations when it comes to acceptance is clearly illustrated 
by the uncorrelated sample, which still appears essentially flat also after the cuts. 

We see, as before, that the 2HDM (11) with large tan/3 differs the most from the V-A 
structure of the SM. Within the statistical variations, it appears difficult to separate between 
the 2HDM (II) with low tan (3 and the SM in any of these channels. 

5 Conclusions 

Top quark spin correlations are predicted by the Standard Model, and this prediction will be 
tested at the LHC. An angular analysis of tt decay products is most efficiently performed in the 
dilepton channel which has a,\ — 1. Unfortunately, the same analysis cannot be carried out in 
the t channel due to the additional neutrinos in the final state which prevents reconstruction of 
the tt CM frame. As we have shown, most of the information about the tt spin correlations can 
however be recovered by looking only in the transverse plane. It would be most interesting to 
investigate the prospects for top spin physics also in the r channel using the tools we suggest. 

Such an analysis would make it possible to extract information about the coupling structure 
of a light charged Higgs boson in the t — > bH + — > br + v T decay mode. The separation in the 
angular variables we discuss between the top decay within the SM, and in the charged Higgs 
decay channel, is hardly of a magnitude interesting for separation of signal and background in 
a discovery phase. Nevertheless, with the good reach to find a light H ± with modest amounts 
of LHC data [18,19], the type of measurements we suggest may still well be considered early 
in comparison to other physics topics beyond the Standard Model. 
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